Exact longitudinal plasmon dispersion relations for one and two dimensional Wigner 

crystals 



O 

o 



(N 



i 

u 

CO 



I 

C 

o 
o 



> 

00 
00 

o 

(N 

^f* 
O 

o 



Shimul Akhanje^B 

Department of Physics, UCLA, Box 951547, Los Angeles, CA 90095-1547 
(Dated: February 1, 2008) 

We derive the exact longitudinal plasmon dispersion relations, u(k) of classical one and two 
dimensional Wigner crystals at T = from the real space equations of motion, of which properly 
accounts for the full unscreened Coulomb interactions. We make use of the polylogarithm function 
in order to evaluate the infinite lattice sums of the electrostatic force constants. From our exact 
results we recover the correct long-wavelength behavior of previous approximate methods. In ID, 
tu{k) ~ jfcj log 1 ' 2 (l/fc), validating the known RPA and bosonization form. In 2D uj(k) ~ \fk. 
agreeing remarkably with the celebrated Ewald summation result. Additionally, we extend this 
analysis to calculate the band structure of tight-binding models of non-interacting electrons with 
arbitrary power law hopping. 

PACS numbers: 



I. INTRODUCTION 

In the realm of one-dimensional(lD) physics, many 
analytical solutions for the interacting electron gas ex- 
ist at various regimes of electron-electron interaction 
strength, both repulsive and attractive, often revealing 
non-trivial fermionic instabilities 1 . Wigner crystals are 
one such phase, existing in the extreme limit of electronic 
correlations 2 . Wigner crystallization is also known to oc- 
cur even in classical systems in higher dimensions for suf- 
ficiently strong unscreened Coulomb repulsion, however 
for the specific case of ID electrons, true long-ranged or- 
der is not possible even with longer ranged interactions. 

Unlike Fermi liquid systems, for Wigner crystals, the 
precise mathematical and experimental behavior of the 
elementary excitations are still far from completely un- 
derstood. The charged collective modes at T = are one 
such entity, where its essential features are captured by 
a classical description^. The spin degrees of freedom also 
lead to another type of collective mode, of which is im- 
portant in describing the magnetic properties. However, 
for simplicity we assume that the spin wave transport is 
sufficiently decoupled from the plasmon propagation. 

The effects of long range interactions are often ne- 
glected in many discussions of Wigner crystals even 
though one would expect in real physical systems a com- 
plete breakdown of screening. So far, the complete plas- 
mon dispersion and more importantly, its precise long 
wavelength behavior has not been deduced analytically. 
We present a simple derivation of the plasmon disper- 
sion relation in a ID Wigner crystal phase that exploits 
the summability properties of power law interactions that 
are unique to ID systems. Morevover, the correct behav- 
ior of the two dimensional(2D) longitudinal eigenmode is 
also calculated exactly by performing the neccessary sin- 
gle summation, of which produces excellent agreement 
with known numerical results. Additionally, as an ex- 
ample of the utility of this form of analysis and owing 
to a similar mathematical structure we also apply these 
methods to calculate the band structure of an electronic 



tight-binding model with power law interactions. This is 
shown in the appendix. 



II. THE DISPERSION RELATIONS 

Let us begin with a classical ID array of L particles 
interacting with unscreened, long-range Coloumb forces, 
described by the Hamiltonian: 



H = 



L 2 

Pi 

^ 2m* 

2—1 c 



(1) 



In the above equation we have defined Pi and m* as the 
particle momentum and effective mass of the ith particle, 
respectively. Our system lives in the low density regime 
in which the Coulomb interactions are much larger than 
the kinetic energy resulting in the crystalline ordering of 
the particles with the spatial coordinates {xi} separated 
with lattice constant a. Moreover, the system is stabi- 
lized by a positive Jellium neutralizing background. The 
eigenvalue equation of interest follows directly from the 
classical equations of motion for the total force acting on 
each particle. This is given by the following expression: 



m P uj 2 u 



(x) + <t>x, x >u{x') = (2) 

x'^x 



where u{x) is the displacement of a lattice site from equi- 
librium and, 



4>x,x' = lj{x- x' f 



(3) 



is the electrostatic force constant between two particles in 
the array, of which is the second derivative of the interac- 
tion potential. As a consequence of the periodic ordering 



2 



and translational symmetry present in the system we can 
assume the eigenfunctions have the form, 



u(na) oc exp[i(kna — ut)] 

of which ft is the Fourier component and n 
We substitute equations and (H]) into @, 



r=l 



sin[fcr/2] 



(4) 

1,2,3.... 

(5) 



where we have defined r = na and have set m e = eo = 1 
for simplicity. Thus, it is our primary task to carry out 
the infinite summation of equation ([5]). Hitherto, the 
most common approaches have been approximate, us- 
ing methods such as the Ewald summation technique 4 , 
where a solution is presented in terms of rapidly converg- 
ing sums. We make use of the Polylogarithm function 
Li n (z) also known as the de Jonquires function, defined 



where Cl n (x) are Clausen functions 5 for a given n. It 
is known from functional analysis that certain Clausen 
functions have an exactly summable representation for 
arguments in a restricted range 5 -. In particular for < 
ft < 2tt, 



„. ., v 7T 2 7rfc k 2 

CUk) = 1 

2W 6 2 4 



(11) 



Apparently, the periodicity of our system guarantees that 
the values of k are restricted to the first Brillouin zone. 
Therefore, a more convenient representation of the 2D 
longitudinal plasma dispersion relation becomes: 



L02D(k) oc 



7T ft 



A' 2 

T 



(12) 



Apparently the long wavelength behavior reduces to: 



Li n (z) 



T- 

k=l 



(6) 



This definition may be extended to all of the complex 
plane through analytic continuation, therefore we apply 
equation ^ to the summation of equation ([5]) yielding, 



&1D 



(fc)«E 



sinfftr] 



r=l 
1 



(7) 



{-Li 3 (e- lk ) - Li 3 (e* k ) + 2((3)) 



where £(x) is the Reimann zeta function. At long wave- 
lengths the polylogarithms can be expanded to the lowest 
order, yielding: 



^(fc^lfcllog^l/lfel] 



(8) 



Let us turn our attention to the two dimensional case, 
for which equation ([2]) can be generalized to a double 
summation over a tensor. The dominant contribution to 
the longitudinal eigenmode is the following sum: 



U2o(k) oc ^ 
r=l 
1 

~ 12 



sinffcr]' 



(n 2 ~ 3Li 2 (e- lk ) - 3Li 2 (e ik )) 



(9) 



We can further simplify this expression by making use of 
the following identity, 

n (*\ - / ¥{ Li n{e~ tx ) - Li n (e lx )} n - even , 
n( ' ~ \ l[Li n (e- lx ) + Li n (e lx )} -» n - odd 1 ' 



U)2D(\k\) OC 



(13) 



If we place this derivation in the context of earlier 
work, this classical result can be compared to the at- 
tempts by other authors using a quantum-mechanical 
treatment of charged collective modes with long-range 
interactions. Until now, no exact analytical results exist 
for the classical plasmon dispersion relations of Wigner 
crystals in any dimension. Gold and Ghazali 6 examined 
a correlated quasi-lD electron system by using the Ran- 
dom Phase Approximation(RPA). In the RPA treatment 
the authors remedy the diverging Fourier transform of 
the \ jr potential in ID by phenomenologically adding a 
small but finite system width d that leads to a logarithmic 
part of the interaction, separating the short-ranged be- 
havior from the long-ranged one. The resulting charged 
modes have the following dispersion, 



oj RPA {k) « ^^Iftlfog 1 / 2 ( M ) 

where vr is the Fermi velocity. Clearly a notable dif- 
ference between our classical result ([7|) and the RPA re- 
sultlfnj is the logarithmic singularity in the limit d — > 0, 
of which is a direct consequence of the authors consider- 
ing a quasi-lD system rather than the purely ID system 
that we have just discussed. Although our exact result 
contains extra dispersive curvature at values of k near the 
Brillouin zone boundary, our classical summation tech- 
nique agrees with the RPA result's long- wavelength be- 
havior. 

Still, others have attempted to describe ID WC behav- 
ior in the limit of an elastic Hamiltonian such as in the 
Luttinger liquid phase(Schulz, 1993) 7 . The Luttinger liq- 
uid arises as an instability in a ID electron system with 
strong short ranged interactions, captured by the Hub- 
bard Hamiltonian. Again they contend that the true long 
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with unscreened Coulomb interactions. Our analysis in- 
troduces mathematical methods for analytically evaluat- 
ing a certain class of lattice sums that have traditionally 
been performed numerically. Recently, Wigner crystals 
in lower dimensions have experienced a resurgence in in- 
terest by both experimentalists and theorists alikeiS, in 
systems such as low density quantum wires and various 
soft condensed matter systems^ 1 . Furthermore, the pre- 
cise wavevector dependence of the eigenfrequencies have 
important consequences for many physically pertinent 
quantities such as the dynamical response functions, and 
thus the results of this paper provide some degree of ana- 
lytical control in future investigations of Wigner Crystals 
that emphasize the role of long-ranged interactions. 



FIG. 1: A comparison of the exact ID classical dispersion 
with the RPA result. 
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FIG. 2: A comparison of the exact 2D longitudinal dispersion 
with the Ewald result. 



range behavior of the Coulomb interaction is not impor- 
tant aside from these minor logarithmic correction factors 
that modify the elastic modes to produce a dispersion re- 
lation that has the same behavior as the RPA result (fH)) . 
A comparison is show in Figure [TJ 

In 2D, until now there have been no analytical results 
for the precise behavior of the longitudinal plasma eigen- 
mode. The Ewald technique has been numerically im- 
plemented by previous authors^. They discovered an 
unusual uj(k) ~ \fk dependence arising strictly from the 
long-ranged interactions, as it is known that short-ranged 
interactions produce a linear phonon dispersive form. We 
compare our closed form result with the Ewald summa- 
tion technique in Figure^ Evidently, there is excellent 
agreement. 



III. CONCLUDING REMARKS 
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APPENDIX A: THE BAND STRUCTURE FOR 
POWER LAW HOPPING 

The tight-binding Hamiltonian of interest takes the fol- 
lowing form, 



ff = £ £ oN)<*l + £MN)01 + b)<*l) ( A1 ) 

where (ij) denotes a sum over all unique pairs of lattice 
sites i and j. For simplicity we have assumed only one 
orbital per site. The hopping matrix elements t%j have 
the power law dependence, 



to 







(A2) 



Evidently, the Hamiltonian satisfies the time- 
independent Schrodinger equation H\ipk) — E(k)\tpk)- 
Furthermore, for translationally invariant systems, we 
can make use of the Fourier transform, having the form 

\tpk) = E elM ' IO: with 2 = 1, 2, ... . We can apply this 

v 

ansatz along with equation (|Aip to the Schrodinger 
equation to yield the following expression, 



io(cos[fcr] — 1) 



= ]T e ifci ' |0 {t {Li p [j k ] + Li [e- ik }) + e } 



(A3) 



We have derived the complete dispersion relation for 
a longitudinal plasmon in ID and 2D Wigner crystals 



where we have defined r = \i — j\ and we have also ab- 
sorbed various constants into the definition of to and eg. 
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Again we made use of the polylogarithm function Li n (z) 
and following similar steps as before our final expression 
for the energy dispersion relation reduces to: 

(k) = t (Lip [e lk ] + Lip [e- ik ] ) + e ( A4) 

Equation (|A4I) is general to a particular exponent /3, and 
one must properly consider the odd and even cases in 



order to determine when it is appropriate use a particular 
Clausen function. More specifically, for odd powers of j3 
one may use the equation (|10p reducing the final energy 
dispersion to: 

E^\k) =2t Clp[k]+e (A5) 
for (3 = 2n + l; n=,0,l,2. ... 



* shimul@physics.ucla.edu 

1 T. Giamarchi, Quantum Physics in One Dimension 
(Clarendon Press, Oxford, 2003). 

2 E. Wigner, Phys. Rev. 46, 1002 (1958). 

3 R. Chitra, T. Giamarchi, and P. L. Doussal, Phys. Rev. B 
65 (2002). 

4 P. P. Ewald, Ann. Phys. 64, 253 (1921). 

5 M. Abramowitz and I. A. Stegun, Handbook of Mathemat- 
ical Functions (Dover Publications, New York, 1970). 

6 A. Gold and A. Ghazali, Phys. Rev. B 41, 7626 (1990). 



7 H. J. Schulz, Phys. Rev. Lett. 71, 1864 (1993). 

8 L. Bonsall and A. A. Maradudin, Phys. Rev. B 15, 1959 
(1977). 

9 R. S. Crandall, Phys. Rev. A 8, 2136 (1973). 

10 P. Glasson, V. Dotsenko, P. Fozooni, M. J. Lea, W. Bailey, 
G. Papageorgiou, S. E. Andresen, and A. Kristensen, Phys. 
Rev. Lett. 87, 176802 (2001). 

11 B. I. Shklovskii, Phys. Rev. Lett. 82, 3268 (1999). 



